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By formally diagonalizing with accuracy h the Hamiltonian of electrons in a crystal subject to 
electromagnetic perturbations, we resolve the debate on the Hamiltonian nature of semiclassical 
equations of motion with Berry-phase corrections, and therefore confirm the validity of the Liouville 
theorem. We show that both the position and momentum operators acquire a Berry-phase depen- 
dence, leading to a non-canonical Hamiltonian dynamics. The equations of motion turn out to be 
identical to the ones previously derived in the context of electron wave-packets dynamics. 

PACS numbers: 



The notion of Berry phase has found many applica- 
tions in several branches of quantum physics, such as 
atomic and molecular physics, optic and gauge theories, 
and more recently, in spintronics, to cite just a few. Most 
studies focused on the geometric phase a wave function 
acquires when a quantum mechanical system has an adi- 
abatic evolution. It is only recently that a possible in- 
fluence of the Berry phase on semiclassical dynamics of 
several physical systems has been investigated. It was 
then shown that Berry phases modify semiclassical dy- 
namics of spinning particles in electric 111 and magnetic 
fields [2J, as well as in semiconductors y|- In the above 
cited examples, a noncommutative geometry, originating 
from the presence of a Berry phase, which turns out to be 
a spin-orbit coupling, underlies the semiclassical dynam- 
ics. Also, spin-orbit contributions to the propagation of 
light have been the focus of several other works 0, 0, la] j 
and have led to a generalization of geometric optics called 
geometric spinoptics Q. 

Semiclassical methods in solid-state physics have also 
played an important role in studying the dynamics of 
electrons to account for the various properties of metals, 
semiconductors , an d insulators Q . In a series of papers 
[1,0 (see also 0), the following new set of semiclassi- 
cal equations with a Berry-phase correction was proposed 
to account for the semiclassical dynamic of electrons in 
magnetic Bloch bands (in the usual one-band approxi- 
mation) 

r = d£(k)/hdk-k x 6(k) 
hk = -eE - er x B(r) (1) 

where E and B are the electric and magnetic fields re- 
spectively and £ (k) = £o(k)— m(k).B is the energy of the 
nth band with a correction due to the orbital magnetic 
moment 0. The correction term to the velocity — k x 
with 0(k) the Berry curvature of electronic Bloch state 
in the nth band is known as the anomalous velocity pre- 



dicted to give rise to a spontaneous Hall conductivity in 
ferromagnets ■ For crystals with broken time-reversal 
symmetry or spatial inversion symmetry, the Berry cur- 
vature is nonzero Eqs^were derived by consider- 
ing a wave packet in a band and using a time-dependent 
variational principle in a Lagrangian formulation. The 
derivation of a semiclassical Hamiltonian was shown to 
lead to difficulties in the presence of Berry-phase terms 
9]. The apparent non-Hamiltonian character of EqsQ] 
led the authors of ^3 to conclude that the naive phase 
space volume is not conserved in the presence of a Berry 
phase, thus violating Liouville's theorem. To remedy this 
situation these authors introduced a modified density of 
state in the phase space D(r, k) = (27r)~ d (l + eB.Q/h) 
such that D(r,k)<irdk remains constant in time. 

This point of view was immediately criticized by sev- 
eral authors. In particular, by relating the semiclassical 
dynamics of Bloch electron to exotic Galilean dynamics 
introduced independently in the context of noncommuta- 
tive quantum mechanics C. Duval et al. found 
that Eqs^&re indeed Hamiltonian in a standard sense, 
restoring the validity of Liouville's theorem when the cor- 
rect symplectic volume form is used. This approach, rely- 
ing on a symplectic structure on a classical Hamiltonian 
formulation, though very elegant, does not stem from the 
quantum Hamiltonian for electrons in a solid and is con- 
sequently not widely known in the solid-state physicists 
community. Additionally, the role of the Berry phase is 
hidden in this approach. In a different but related work 
|l5|. the Hamiltonian nature of semiclassical equations 
of motion of Dirac electrons in electromagnetic field with 
Berry-phase corrections (in this case it corresponds to a 
spin-orbit coupling) was established. 

This letter presents an alternative approach for the 
derivation of the equations of motion of an electron in 
magnetic Bloch bands, based on a direct semiclassical di- 
agonalization of the full quantum Hamiltonian. We show 
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that both dynamical variables (r. k) in Eqs. are not 
the usual Galilean operators, but new covariant opera- 
tors denned in a particular nth Bloch band, and includ- 
ing Berry gauge potentials. These potentials induce non- 
canonical commutation relations between the covariant 
variables. In our context, the equations of motion are 
given by the tandard dynamical laws hdv/dt = i[r,H] 
and dk/dt = i[k, H] leading to Eqs^in a semiclassical 
approximation. Our approach thus reveals the Hamilto- 
nian nature of Eqs^and confirms the importance of the 
Berry phase on the semiclassical dynamics of Bloch elec- 
trons. The origin of the density of state D(r, k) is then 
obvious; it is simply equal to the Jacobian of the transfor- 
mation between the canonical variables (R, K) and the 
covariant ones (r, k) , as already found in the context of 
the Dirac equation in (lif. 

It should be noted that focusing on a Hamiltonian for- 
malism for electrons in solids in order to account for 
the anomalous velocity was first initiated by Adams and 
Blount 0|, who showed that this term arises from the 
noncommutability between the components of the intra- 
band position operator, which acquires a Berry-phase 
contribution. But their approach does not lead to the 
correct EqsQ for electrons in magnetic Bloch waves, as 
they missed the Berry phase dependence of the intraband 
momentum operator. A similar Hamiltonian approach 
has also been realized for arbitrary spinning (massive and 
massless) particles in an electric field |l| and extended to 
the case of Dirac electron in an arbitrary electromagnetic 
field HE1- 

The common feature of these Hamiltonian 
formulations is that a noncommutative geometry under- 
lies the algebraic structure of both coordinates and mo- 
menta. Actually, a Berry-phase contribution to the coor- 
dinate operators stems from the representation where the 
kinetic energy is diagonal (Foldy-Wouthuysen or Bloch 
representation). The components of the coordinate be- 
come noncommutative when interband transitions are ne- 
glected (adiabatic motion). 

Consider an electron in an crystal lattice perturbated 
by the presence of an external electromagnetic field. As is 
usual, we express the total magnetic field as the sum of a 
constant field B and small nonuniform part <5B(R). The 
Hamiltonian can be written H = Hq — e0(R), with Hq 
the magnetic contribution (0 being the electric potential) 
which reads 

H = ( ^- + eA(R) + eSA(R)) +V(R) (2) 
\ Zm J 

where A(R) and (5A(R) are the vectors potential of 
the homogeneous and inhomogeneous magnetic field, re- 
spectively, and V(R) the periodic potential. The large 
constant part B is chosen such that the magnetic flux 
through a unit cell is a rational fraction of the flux quan- 
tum h/e. The advantage of such a decomposition is 
that for (5A(R) = the magnetic translation operators 
T(Ri) = exp(iK.Ri), with K the generator of transla- 



tion, are commuting quantities allowing to exactly diag- 
onalize the Hamiltonian and to treat <5A(R) as a small 
perturbation. The state space of the Bloch electron is 
spanned by the basis vector |n, k >= |k > ®|n > with 
n corresponding to a band indice. In this representa- 
tion K|ri, k >= k|n,k > and the position operator is 
R = id/dH., which implies the canonical commutation 
relation [Ri, Kj] = idij. 

We first perform the diagonalization of the Hamilto- 
nian in Eq|5| for <5A ~ by an unitary matrix U (K) 
(whose precise expression is not necessary for the deriva- 
tion of the equations of motion) such that UHU + = 
£(K) — e0(f/RJ7 + ), where £(K) is the diagonal energy 
matrix made of elements £ n (K) with n the band indicc. 
Whereas the quasi-momentum is invariant through the 
action of U, e.g., k = U"KU + = K, the position operator 
becomes: 

v = umj + = n + iUchjj + (3) 

in the new representation owing to the fact that 
[Ri,Kj\ = iSij. In the adiabatic or one-band approx- 
imation, which consists of neglecting interband transi- 
tions, one has to project the position coordinate (the 
momenta operator is diagonal and so invariant by con- 
struction) on a certain band such that the nth intra- 
band position operator r„ = V„ (r) reads r„ = R + A n . 
The quantity A n = i r P n (Ud] s _U + ) is a Berry connection, 
as it can be readily shown that its matrix elements are 
given by A n {k) = i < u„(k)|<9k|w„(k) >, where we used 
U + (k)\n >= |u n (k) > with |u n (k) > the periodic part 
of the magnetic Bloch waves. The price to pay when 
considering the one-band approximation is that the alge- 
bra of the coordinates becomes noncommutative (as we 
consider only one band, we drop the index n) 

[r\r j ] = i6 y '(k) (4) 

with 6 iJ (k) = aU 3 (k) - c*U l (k) the Berry curvature. 
Observe that the replacement of k by p/h shows that 
O lJ (p) is actually of order h 2 . In the one-band approx- 
imation the full Hamiltonian, including the electric po- 
tential, is now given by 

V n {UHU+) = £(k) - e^(r). (5) 

Due to the Berry connection in the definition of the posi- 
tion operator, the equations of motion should be changed. 
But to compute commutators like [r k , </>(r)] , one re- 
sorts to the semiclassical approximation [r fc ,(/>(r)] = 
idi<j)(r)Q kl + 0(h) leading to the following semiclassical 
equations of motion 

r = d£{k)/hdk -kx 9(k), hk = -eE (6) 

where E is the external electric field. Whereas the mo- 
mentum equation of motion is the usual one, the velocity 
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operator acquires an anomalous contribution due to pres- 
ence of the Berry curvature. Notice that the contribution 
of the magnetic field stems only from the presence of the 
Berry curvature through the band structure. This equa- 
tion was first derived by Adams and Blount using a 
similar approach, and later rederived by Niu and cowork- 
ers [1, 13 by looking at the dynamics of wave packets 
from a Lagrangian formalism. In the following, we will 
extend our approach to carry out a semiclassical diago- 
nalization of the full electromagnetic Hamiltonian (with 
<5A(R) 7^ 0). Contrary to the work of [16J, we show 
that the momentum also acquires a Berry-phase contri- 
bution leading to different semiclassical equations of mo- 
tion. These last ones turn out to be those derived first 
m 0E} (also Duval et al. Q in another context). Our 
rigorous approach has the merit to show without ambigu- 
ities that the equations of motion are indeed Hamiltonian 
in the standard sense. 

The diagonalization of the Hamiltonian in the presence 
of an arbitrary magnetic field is now the focus of the rest 
of the paper. Consider first the Hamiltonian Eq|21in the 
absence of an electric field and set K = K + e<5A(R)/fi,. 
As the flux (5B on a plaquette is not a rational multiple of 
the flux quantum, we cannot diagonalize simultaneously 
its components Ki since they do not commute anymore. 
Actually 



obtain 



k = UKU+ = K + Ud &j U^ 



h[k\K J ] = -iee ijk SB k (TL) 



(7) 



As a consequence of this non-commutativity, we just 
aim at quasi-diagonalizing our Hamiltonian at the semi- 
classical order (with accuracy K). To perform this ap- 
proximate diagonalization U HU + with accuracy H we 
first consider the limiting case of a constant potential 
<5A(R) = <5Ao (this is obviously a formal consideration). 
Clearly, the Hamiltonian in Eq(21 is diagonalized by the 
matrix U (SA) = U(K + eSA/h), as we have just shifted 
the momentum K. To diagonalize Eq. [3 in the general 
case, we now consider a unitary matrix U (K+e<5A(R) / h) 
which has the same series expansion as U (8 A (R)) when 
R is considered as a parameter commuting with K. Of 
course, this matrix is not unique, due to the noncom- 
mutativity of K and R, but it can be shown that the 
various choices lead to the same projected Hamiltonian. 
Note that in the sequel, a small 5 A perturbation, which 
preserves the band structure determined previously is as- 
sumed, i.e. < n\5A\m > = for m n. Before im- 
plementing effectively the canonical transformation on 
the Hamiltonian, it appears more convenient to imple- 
ment first the canonical transformation on the dynami- 
cal operators. Therefore, in the new representation the 
position operator is again given by r = R + iUd^U + . 
As before, the projection on a band defined the nth 
intraband position operator r„ = R + „4„(K), with 
A n (K) = V n {Ud-g L U + ) a new Berry connection. 

The pseudo-momentum K is no more invariant as we 



= K — ieUd il U + x SB(R)/h (8) 

The nth intraband momentum operator k n = P n (k) is 
then 

HK = HK - eAn(K) x £B(R) (9) 

which at the order h can also be written 

hk n ~ HK - &A(k n ) x SB(r n ) + 0{h 2 ) (10) 

This new contribution to the momentum has been over- 
looked before in the work of Adams and Blount 0] but is 
crucial for the correct determination of the semiclassical 
equations of motion of an electron in a magnetic Bloch 
band. 

The commutation relations between the components of 
the intraband momenta are therefore given by (at leading 
order) 

h = -iee tlk 8B k {r n )+ie 2 e tpk 8B k e ]ql 5BiQ pq /h 

(11) 

with 9^(k„) = <9M J (k„) - <9M 4 (k„) the Berry curva- 
ture. The commutation relation between position and 
momentum can be computed leading to 

<, Hki\ = ih5 ij + iee ]lk SB k (r n )e u (k n ) (12) 

The third useful commutator is as in Eq0] given by 

=iB(Kr (13) 

at leading order. The set of nontrivial commutations re- 
lations given by Eqs^I 1^1 QUI is the same as the one 
deduced in [l5| in the context of the Dirac electron us- 
ing an approximate explicit Foldy Wouthuysen transfor- 
mation. This shows that a common structure underlies 
the quasi-diagonalization of general quantum Hamiltoni- 
ans in the presence of electromagnetic fields ll7|. In the 
present case, the approximate diagonalization UHU + is 
performed by formally expanding U and H in a series of 
K and <5A(R). The recombination of the series includes 
corrections of order H due to the noncommutativity of K 
and R. In doing so, we arrive at the following expression 



UHU+ = £ { k )--H 
ie 



£(K),A(K)1 E ijk SB k (K)Aj(K) 



— A,-(K) [£(K),MK)\ e^ k SB k (K) 
which after projection on the nth band can be written: 
V n (UHU+) = £ n (k„) - M(K).6B(r n ) + 0(h 2 ) (14) 
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with M(K) = £>„(§f £(K),A(K) x A{K)) the magne- 
tization. This term can also be written under the usual 
form in the (K,n) representation (lif : 



Mi 



2h 



(15) 



We mention that this magnetization (the orbital mag- 
netic moment of Bloch electrons) , has been obtained pre- 
viously in the context of electron wave packets dynamics 



Notice that because a semiclassical computation was 
considered here, we kept only terms of order h. As SA is 
small, we chose to neglect terms of order MA 2 . But, as 
we do not consider a perturbation expansion, our method 
keeps all contributions of order 5A n . In a perturbation 
expansion, instead of evaluating U HU + , one would com- 
pute U(K)HU(K)+ = £(K) + USHU+ (and neglect all 
terms of order higher than <5 A ). In this representation 
the position operator is still given by Eq[3] but K is in- 
variant. But doing so would lead us to neglect contri- 
butions of order h that are fundamental for the correct 
determination of the equations of motion. A perturba- 
tion expansion is then not allowed here. 

The commutation relations EqsEl 113 EUl together 
with the semiclassical Hamiltonian of the Bloch elec- 
tron in the full electromagnetic field E n (k n ) — (f>(r n ) with 
E n (k n ) = £ n (k n ) — A / ((k,i).5B(r„), allow us to deduce 
the semiclassical equations of motion. Dropping now the 
index n we have: 

r = dE(k)/hdk-k x 9(k) 
hk = -eF,-erxSB(r)-M.dSB/dr (16) 

These equations differ from the ones derived in [la. but 
are exactly the same equations of motion as in apart 
from the magnetization contribution (which should also 
be present in 0). It is also clear that the correct volume 
form in the phase space (r, k) has to include the Jaco- 
bian D(r, k) = (l + eJB.O/fi.) of the transformation from 
(R, K) to (r,k). This Jacobian is precisely the density 
of state introduced in ^3] , m order to ensure the validity 
of the Liouville theorem. As a consequence, and by com- 
paring Easfilandllfilwe can conclude that the operators 
(r, k) correspond to the dynamical variables in Eqs. ^ 
denoted x c and q c in Ref. [9|. The variable x c is the 
wave-packet center, and q c the mean wave vector. The 
equations of motion for x c and q c were obtained, using 
a time-dependent variational principle in a Lagrangian 
formulation . It was then found that the derivation of 
a semiclassical Hamiltonian presents some difficulties in 
the presence of Berry-phase terms. Actually, as explained 
in |9J] , this derivation requires the knowledge of the com- 
mutation relations between x c and q c (a re-quantization 



procedure) , but these relations cannot be found from the 
Lagrangian formulation. One of the advantages of our 
approach is to show that these commutation relations 
are in fact a direct consequence of the semiclassical diag- 
onalization of the quantum Hamiltonian. 

In summary, our semiclassical diagonalization of the 
electromagnetic Bloch Hamiltonian leads to a well de- 
fined semiclassical Hamiltonian with Berry-phase correc- 
tions. The resulting semiclassical equations turn out to 
be the ones obtained previously from a semiclassical La- 
grangian formalism 9] . When the correct dynamical vari- 
ables are used the Liouville theorem is restored. More- 
over, the present approach also confirms the result of 
Duval et al. and Bliokh Jlgj about the Hamiltonian 
nature of these semiclassical equations of motion with 
Berry-phase corrections, which is a hotly debated sub- 
ject. 

We would like to thank Aileen Lotz for a critical read- 
ing of the manuscript, and one referee whose pertinent 
questions allowed us to improve the present article. 
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